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ABSTRACT 

We present first results from a 325 ks observation of the Seyfert 1 galaxy MCG-6-30-15 with 
XMM-Newton and BeppoSAX. The strong, broad, skewed iron line is clearly detected and is 
well characterised by a steep emissivity profile within 6r g (i.e. QGM/c 2 ) and a flatter profile 
beyond. The inner radius of the emission appears to lie at about 2r g , consistent with results 
reported from both an earlier XMM-Newton observation of MCG-6-30-15 by Wilms et al. 
and part of an ASCA observation by Iwasawa et al. when the source was in a lower flux state. 
The radius and steep emissivity profile do depend however on an assumed incident power-law 
continuum and a lack of complex absorption above 2.5 keV. The blue wing of the line profile 
is indented, either by absorption at about 6.7 keV or by a hydrogenic iron emission line. The 
broad iron line flux does not follow the continuum variations in a simple manner. 
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1 INTRODUCTION 

The Seyfert 1 galaxy MCG-6-30-15 has played an important role 
in studies of accretion onto black holes due to the presence of a 
broad, skewed iron line in its X-ray spectrum (Tanaka et al. 1995). 
The shape of the line seen with ASCA is consistent with emission 
from the surface of an accretion disc extending from about 6 to 
more than 40 gravitational radii (6 — 40r g ; r g = GM/c 2 ) inclined 
at about 30 deg to the line of sight (Fabian et al. 1995). Occasionally 
the red (i.e. lower energy) wing of the line is seen to extend below 
4 keV (Iwasawa et al. 1996; 1999). This can be explained by the 
disc extending within 6r g which may imply the black hole must be 
rapidly spinning. The presence of the broad iron line in MCG-6- 
30-15 has been confirmed by BeppoSAX (Guainazzi et al. 1999), 
XMM-Newton (Wilms et al. 2001) and Chandra (Lee et al. 2002). 

Here we present preliminary results from a long 325 ks obser- 
vation of MCG-6-30-15 made with XMM-Newton. The source was 
at a similar flux level to the previous ASCA observations, and about 
70 per cent brighter than during the earlier 100 ks XMM-Newton 
observation reported by Branduardi-Raymont et al. (2001) and 
Wilms et al. (2001). Simultaneous observations were made with 
BeppoSAX, providing coverage from ~0. 2-100 keV. The present 
work focuses on the spectrum above 2.5 keV and the iron-K line 
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features; absorption and emission features below 2 keV due to oxy- 
gen and other elements will be discussed more fully in later work. 



2 DATA REDUCTION 

MCG-6-30-15 was observed by XMM-Newton (lansen et al. 2001) 
over the period 2001 luly 31 - 2001 August 5 (rev. 301, 302 and 
303), during which all instruments were operating nominally. The 
present analysis is restricted to the data from the European Photon 
Imaging Cameras (EPIC). Both the EPIC MOS cameras (Turner et 
al. 2001) and the EPIC pn camera (Striider et al. 2001) were oper- 
ated in small window mode and used the medium filter. Extraction 
of science products from the Observation Data Files (ODFs) fol- 
lowed standard procedures using the XMM-Newton Science Anal- 
ysis System v5.2 (SAS). 

The EPIC data were processed using the standard SAS pro- 
cessing chains. Source data were extracted from circular regions of 
radius 30 arcsec from the processed MOS and pn images and back- 
ground events were extracted from regions in the small window 
least effected by source photons. These showed the background to 
be relatively low and stable throughout the observation, with the ex- 
ception of the final few ks of each revolution where the background 
rate increased. Data from these periods were ignored. The total 
amount of "good" exposure time selected was 315 ks and 227 ks 
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Figure 1. Combined MOS spectrum (the data were combined for plotting 
purposes only) and ratio of data to a power-law model joining the 2-3 keV 
data and 7.5-10 keV data. As this is not a realistic model for the continuum 
the residuals should be considered merely as representative of the spectral 
complexity. The inset panel shows a close-up of the iron line region with 
the 6.9 keV feature marked. 



(specifically line and edge energies) are quoted for the rest frame 
of the source. 

The EPIC MOS data show a spectral form remarkably sim- 
ilar to those seen in earlier ASCA observations, as illustrated in 
Fig |l| (cf. Fig. 1 of Iwasawa et al. 1996). The residuals around 
6 keV have been interpreted as a broadened and redshifted iron 
Ka emission from a relativistic accretion disc (Tanaka et al. 1995). 
Only the XMM-Newton EPIC spectrum above 2.5 keV is exam- 
ined here, to limit the effects of the complex soft X-ray absorption 
(see also Lee et al. 2001), a detailed treatment of which is beyond 
the scope of this Letter. This energy range also excluded the re- 
gion around the detector Si K-edge (1.84 keV) and the mirror Au 
M-edge (2.3 keV). To better constrain the continuum, and in par- 
ticular the Compton reflection component (e.g. George & Fabian 
1991; Nandra & Pounds 1994), the BeppoSAX MECS and PDS data 
were also fitted simultaneously. The normalisations of the MOS 
and MECS data were allowed to vary independently, to account for 
differences in the absolute calibration, with the PDS normalisation 
tied to 0.86 times the MECS value (Fiore et al. 1999). Galactic ab- 
sorption (A^h = 4.06 x 10 20 cm" 2 ; Dickey & Lockman 1990) was 
included in all the models discussed below. 



for MOS and pn, respectively. (The lower pn exposure is due to the 
lower "live time" of the pn camera in small-window mode, ~71 per 
cent). 

The ratios of event patterns as a function of energy showed 
there is negligible pile-up in the pn data but the MOS data suffer 
slightly from pile-up. A comparison of the data extracted includ- 
ing and excluding events from the central 10 arcsec of the source 
region showed only a slight flattening (AF ~ 0.06) in the piled- 
up spectrum. This small effect was ignored in order to increase the 
number of source counts. 

Events corresponding to patterns 0-12 (single-quadruple 
pixel) were extracted from the MOS cameras and patterns 0-4 (sin- 
gle+double) were used for the pn spectral analysis, after checking 
for consistency with the pattern (single pixel event) data. The total 
number of photons extracted from the source region was 6.23 x 10 6 
for the pn and 2.35 x 10 6 for each MOS camera. The BeppoSAX 
data were processed in the standard manner (Fiore et al. 1999). 



3 SPECTRAL ANALYSIS 

The spectral calibration of the EPIC cameras is still evolving. Using 
the publicly available software and response matrices, the remain- 
ing calibration-dependent residuals over the 2.5-10 keV range con- 
sidered here should be at the <J 5 per cent level (see Ferrando et al. 
2002). However, there remain problems with the model for charge 
transfer inefficiency (CTI) for pn small-window mode data. In or- 
der to mitigate the effects of this, a conservative approach was used 
for the present analysis, with detailed spectral fitting results based 
on only the MOS data. 

The source spectra were grouped such that each bin 
contains at least 20 counts and were fitted using XSPEC 
Vll.l (Arnaud 1996). The spectral responses used were 
rtil_medv9ql9t5r5_all_15 . rsp for MOS1 (and similarly for 
MOS2) and epn_sw2 0_sdY9jnedium. rmf for the pn. The 
quoted errors on the derived model parameters correspond to a 
90 per cent confidence level for one interesting parameter (i.e. a 
A\ 2 = 2.7 criterion), unless otherwise stated, and fit parameters 



3.1 Phenomenological model 

The iron line residuals were initially parameterised using a simple 
model comprising a power law plus Compton reflection, narrow 
(a = 10 eV) Gaussian at 6.4 keV and a broad Gaussian whose en- 
ergy and width were free parameters. The Compton reflector was 
assumed to be a neutral, static slab subtending a solid angle of 27r sr 
at the X-ray source, inclined at i = 30 deg, and having solar el- 
emental abundances (Magdziarz & Zdziarski 1995). The primary 
continuum was a power law with an exponential cutoff at a fixed 
energy of 400 keV (constraints on the reflection and continuum pa- 
rameters are discussed later). This model gave a good overall fit 
with x 2 = 984.1 for 1034 de grees of freedom (dof). The broad 
line component had an energy E = 4.8 ± 0.3 keV, a width of 
a = 1.8 ± 0.2 keV and equivalent width EW = 527 ± 70 eV. The 
narrow, 6.4 keV "core" had EW = 38 ± 10 eV and peaked very 
close to 6.4 keV, with E = 6.40 ± 0.04 keV if allowed to be free. 

Very narrow components to the iron line, unresolvable by 
EPIC, have been observed in the Chandra HETGS spectra of sev- 
eral Seyferts (Yaqoob et al. 2002). They may be identified with 
emission from, e.g. the BLR or molecular torus. Such narrow lines 
are not ubiquitous, however, and allowing the 6.4 keV component 
to have a finite width is resolved with a width corresponding to 
FWHM as 30, 000 km s _1 and EW = 120 eV. The fit is excel- 
lent, the parameters are given in Table [l] (model 1). A preliminary 
analysis of the pn data gave generally consistent results and in par- 
ticular the core of the line was again resolved. See also Lee et al. 
(2002). The large width of the line core suggests that this also arises 
in the disc, rather than at a greater distance (cf. NGC 5548; Yaqoob 
et al. 2001). 



3.2 Relativistic line profiles 

To provide a more physically realistic description of the data the 
line profile of Laor (1991) was used to describe a broad line from 
an accretion disc surrounding a rotating black hole. The disc was 
assumed to be in a low state of ionisation, with the rest energy of the 
line fixed at 6.4 keV. The fits were relatively insensitive to the outer 
radius of the disc, which was therefore fixed at R ou t ~ 400r g . The 
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Table 1. Results of simultaneous fits to the XMM-Newton MOS and BeppoSAX data. Line energies (E) and widths (a) are given in units of keV and equivalent 
widths (EW) are given in eV. Disc radii (e.g. R ln ) are given in units of r g and inclination (i) is given in degrees. The models are described in section pi 
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Continuum 
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i 

38.3 ± 1.2 


Broad Ka 
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554 ± 56 
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EW 
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X 2 /dof 
885.1/1031 
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27.8 ± 1.3 


2.0 ±0.2 




4.8 ±0.7 




41 Q+ 73 


21 ±5 


856.2/1029 



1 This model includes R ICT \ and -Efoid as free parameters, Ap c = 3 and emission from Fe K/3 (see section 




free parameters in the model were: the inner radius of the disc, R ln , 
the inclination, i, and the emissivity as a function of radius, which 
was parameterised as a power law of index q (i.e. R~ q ). Both the 
line and the reflection spectrum were blurred using the kernel from 
the LAOR model. The fit parameters are given in Table [I] (model 2), 
the x 2 is comparable to the broad plus narrow Gaussian model, but 
falls short of the best fitting vj 2 with two broad Gaussians (model 
1). The reason for this is that, with a simple power law emissiv- 
ity function, the Kerr line profile cannot simultaneously model the 
broad, red wing of the line and the narrower core around 6.4 keV. 
The preliminary pn analysis gave comparable results but suggested 
an even stronger red wing to the line. 

3.3 Additional iron K components 

Before exploring the red wing further it should be noted that there 
is some additional complexity in the line region at ~ 7 keV (see 
inset of Fig [j]). Including an iron K-edge in model 2 provided a sig- 
nificant improvement in the fit (x 2 = 908.2/1032), with an edge 
energy of E = 7.38 ± 0.09 keV and a depth of r = 0.11 ± 0.02, 
similar to that discussed by Pounds & Reeves (2002). However, 
such an edge must have a physical origin; if it is due to neutral or 
near-neutral iron along the line of sight, then strong iron-L absorp- 
tion is expected which is inconsistent with the soft X-ray spectrum 
(Lee et al. 2001). If the material is out of the line of sight and the 
edge is due to reflection off distant material, the reflection com- 
ponent must be very strong (R ~ 3), which is then inconsistent 
with the relatively weak, narrow 6.4 keV iron line observed in the 
spectrum. Another possibility is that the edge arises from thick ma- 
terial partially covering the line of sight. This gives a comparable 
fit (x 2 = 894.6/1032) to the edge, with a covering fraction of 
~ 0.16. Using this model the broad line parameters differed only 
slightly from those of model 2, and the broad, red wing of the line, 
extending within 6r g , was still required in the fit. 

Sako et al. (2002) suggested there may be significant ab- 
sorption in the 6.4-6.7 keV range due to inner-shell transitions of 
ionised iron. This possibility was explored by including an absorp- 
tion line in model 2. This gave a good fit (x 2 = 923.9/1031 dof) 
with absorption line parameters E = 6.74 ± 0.05 keV and EW = 
138 ± 35 eV. The relativistic line parameters differ only slightly, 
with the emissivity index changing tog = 6.6± 1.1. 

An alternative interpretation is that the complexity around 
7.0 keV is in an emission component, identified with recombina- 
tion emission by H-like iron. Adding a Gaussian to model 2 at 
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Figure 2. Model including relativistic iron line emission (model 4) and fit 
residuals for XMM-Newton MOS (black) and BeppoSAX MECS (green) and 
PDS (blue). 

E ta 6.9 keV again improved the fit ( X 2 = 956.1/1032 dof) 
and gives best fitting parameters of E = 6.91 ± 0.03 keV and 
EW = 18 ± 6 eV. If a line of this strength is due to gas pho- 
toionised by the observed continuum then it must lie within a radius 
of ~ 2 x 10 16 cm. Since it is not possible to unambiguously deter- 
mine between absorption and emission as the source of complexity 
above 6.5 keV, both are considered in the following subsections. 

3.4 More complex models 

A steep emissivity function was required to model the red wing of 
the line (see also Wilms et al. 2001), but this gave too little flux 
from the outer parts of the disc to model the core of the line. This 
raises the possibility that the emissivity function is not a simple 
power-law. This was explored by using a broken power-law for the 
emissivity function, with emissivity index q ln from R[ n to Rhr and 
index g out from Rhr to R ou t- The fit parameters for this model are 
given in Table]!] (model 3). 

This fit is not physically self-consistent for at least two rea- 
sons. Firstly it ignores the emission from Fe K/3, which is expected 
at 7.05 keV. Secondly, the high EW of the line is not consistent 
with our assumptions about the reflection component, namely that 
it has .Rrcfl = 1, i = 30 deg, solar abundances and is neutral. In 
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Figure 3. Relativistic iron line profile, shown in (-Fi/) flux units, obtained 
from the ratio of the MOS data to best-fitting underlying continuum model 
(model 4), multiplied by the continuum model in flux units (as opposed to 
an "unfolded" plot). The crosses mark the data points and the solid line 
marks the line model. 



such circumstances, considering the observed value of V, the pre- 
dicted EW is only 140 eV (George & Fabian 1991). Two plausi- 
ble possibilities are that the iron abundance in the accretion disc 
is higher than the solar value (see e.g. Lee et al. 1999), or that 
the reflection component is enhanced for some reason, for exam- 
ple due to anisotropy (Ghisellini et al. 1991) or geometry (Fabian 
et al. 2002). 

The high abundance case was examined by fixing the iron 
abundance of the reflector to 3 x solar in model 3. This in fact 
improved the overall fit compared to the solar abundance case 
(X 2 = 877.2/1031 dof). Changes in the disc line parameters 
were negligible. The value of the cutoff energy can affect both 
the high energy spectral shape and the apparent strength of the 
reflection component. Allowing E cut to be free in this fit gave 
E cut = 241±86 8 keV ( X 2 = 875.4/1030 dof). Allowing i? rofl 
to be a free parameter further improved the fit, therefore this was 
also left free in model 4. 

Model 4 comprised a power-law plus cold reflection (with 
overabundant iron, Af c — 3), emission from iron Kq (6.40 keV) 
and K/3 (7.05 keV). The reflection plus line spectrum was then 
blurred using the LAOR kernel (with a broken power-law emissivity 
function). A narrow emission line at E = 6.9 keV (fixed) was also 
included in the model. This provided the best fitting model, with 
a cut-off energy E cut = 122±?£ keV and i? rcfl = 2.2±J;^, the 
other parameters are given in Table |l|. The model and fit residuals 
are shown in Fig. ^ and the "fluxed" line profile based on this fit is 
shown in Fig. pi 

Including an absorption line instead of a 6.9 keV emission 
line gave a slightly poorer fit (\ 2 ~ 863.0/1027 dof) with line 
parameters E = 6.74 ± 0.02 keV, EW = 31 ± 9 eV (consistent 
with the prediction of Sako et al. 2002), and slight changes to the 
relativistic line parameters (in particular 7?br = 3.9 ± 0.7r g ). 



4 SPECTRAL VARIABILITY 

As a first test for spectral variability around the iron line region pn 
spectra were extracted from 10 ks intervals around the minimum 
and maximum flux periods during revolution 301 (see Fig. Q). The 
average continuum fluxes during these intervals differ by > x2. 
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Figure 4. EPIC pn light curve (0.2-10 keV) for revolution 301 showing the 
intervals of low and high flux discussed in section W. 
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Figure 5. Comparison of the high-flux and low-flux pn spectra extracted 
from intervals marked on Fig. W. Top panel: ratio of high-flux to low-flux 
spectrum. Bottom panel: comparison of the difference spectrum (high-low) 
to a power-law model. 



The ratio of the two spectra is shown in the top panel of Fig. ^| and 
shows clear signs of spectral variability. There is a slight depression 
around 6.4 keV, indicating the iron line core is stronger (compared 
to the continuum) in the low-flux spectrum. The bottom panel of 
Fig. |^ shows the residuals from fitting a power-law model to the 
difference spectrum (i.e. high-flux minus low-flux spectrum). The 
difference spectrum is consistent with a power-law in the range 2- 
10 keV and indicates the iron line flux changed little between the 
two intervals. A detailed analysis of the spectral variability prop- 
erties will be given in a forthcoming paper but these first results 
are consistent with the simple two component continuum model 
described in Shih et al. (2002). 

The suggestion that the iron line was not as variable as the 
continuum was borne out by an examination of the RMS spectrum. 
Figure |fj shows the normalised RMS variability spectrum derived 
from the revolution 302 pn data. The RMS spectrum clearly shows 
that the fractional variability amplitude is suppressed at energies 
close to the iron line, compared to surrounding continuum bands 
(see also Inoue & Matsumoto 2001). 

A further point to note is that the RMS spectrum (and the ratio 
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Future work to determine the spin must include emission 
from the immediate plunge region inside the innermost stable or- 
bit (Reynolds & Begelman 1997; Agol & Krolik 2000) as well as 
returning radiation (Martocchia, Matt & Karas 2002). 
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Figure 6. The RMS spectrum based on rev. 302 EPIC pn data. The errors are 
calculated as in Edelson et al. (2002) and should be considered conservative 
estimates of the true uncertainty. 



of high/low spectra) are smooth in the 0.2-1.0 keV range. In par- 
ticular, the absence of any feature at 0.7 keV in the spectral ratio is 
expected if the large drop at that energy (Fig |l|) is due to absorption 
(with constant optical depth). For the drop to be the blue wing of an 
emission line (Branduardi-Raymont et al. 2001), the line intensity 
must be responding linearly to the continuum flux (unlike the iron 
line). An examination of the RGS spectra confirms this result. 



5 DISCUSSION 

A long observation with XMM-Newton and BeppoSAX of MCG-6- 
30-15 in its typical state has again confirmed the presence of the 
broad, skewed iron line. All the models considered above are for- 
mally acceptable (xt < 1.0), and all include emission within 6r g , 
consistent with emission from a disc around a spinning black hole. 

Model 4 provides the best (in a x 2 " sense ) an d most physi- 
cally self-consistent explanation of the data. In this model, the disc 
emissivity is described by a broken power-law in radius, and it is in- 
teresting to note that the break radius occurs at ~ 6r g . Beyond this 
radius the disc has an emissivity profile q out ~ 2.5 and produces an 
iron line with an equivalent width ~ 200 eV (the 5.5 — 6.5 keV core 
of the line shown in Fig. ^), both close to the values expected from 
standard accretion disc models. Within 6r g the emissivity steepens, 
producing the strong low-energy tail to the line emission, also with 
an equivalent width ~ 200 eV, suggesting additional physics within 
this region beyond that expected from standard accretion disc mod- 
els (see Wilms et al. 2001). The strong iron line is consistent with 
an enhanced reflection spectrum and an overabundance of iron or 
ionisation of the disc surface. This last possibility will be examined 
in a later paper. 

The iron line parameters do, however, depend slightly on the 
assumptions made about the complex absorption. In particular all 
the above models formally assume the effects of ionised absorption 
are negligible above 2.5 keV, which is correct for absorption due to 
low-Z elements (e.g. O) but may not be true if there is absorption 
by higher-Z elements (e.g. Si, S). A preliminary analysis allowing 
for the possibility of Si and S edges did not change the requirement 
for a strong red wing to the iron line. A detailed analysis of the RGS 
data will yield constraints on the warm absorption (particularly on 
the low Z elements) and help remove existing degeneracies. 
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